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We show that the controversies on the gauge dependence and the infrared singularity emerged in 

G\ • 

Q\ ■ the generalized factorization approach for nonleptonic heavy meson decays within the framework 

of the operator product expansion can be resolved by perturbative QCD factorization theorem. 
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Abstract 



Gauge invariance of the decay amplitude is maintained under radiative corrections by assuming on- 
shell external quarks. For on-shell external quarks, infrared poles in radiative corrections have to be 
extracted using the dimensional regularization. These poles, signifying nonperturbative dynamics of 
a decay process, are absorbed into bound-state wave functions. Various large logarithms produced 
^ ■ in radiative corrections are summed to all orders into the Wilson and Sudakov evolution factors. The 

remaining finite part gives a hard subamplitude. A decay rate is then factorized into a convolution 
of the hard subamplitude, the Wilson coefficient, and the Sudakov factor with the bound-state 
wave functions, all of which are well-defined and gauge invariant. 



I. INTRODUCTION 

The effective Hamiltonian is the standard starting point for describing the nonleptonic 
weak decays of hadrons. Consider the decay B — ► D + tt~ as an example. The relevant 
effective AB = 1 weak Hamiltonian is 

(jc F 

V2 



ftcs - —7= V cb V* d 



Ci(/X)0i(//) +C2(//)0 2 (/x) 



where 



Oi = {cb) v _ A (du) v _ A , 2 = (db) v _ A (cu) v _ A , (2) 



with {q\q2) v±A = 5i7/u(l ± 75)^2- in order to ensure the renormalization-scale and -scheme 
independence for the physical amplitude, the matrix elements of 4-quark operators have 
to be evaluated in the same renormalization scheme as that for Wilson coefficients and 
renormalized at the same scale \i. 

Although the hadronic matrix element (0(h)) can be directly calculated in the lattice 
framework, it is conventionally evaluated under the factorization hypothesis so that (O(fi)) 
is factorized into the product of two matrix elements of single currents, governed by decay 
constants and form factors. In spite of its tremendous simplicity, the naive factorization 
approach encounters two principal difficulties. First, it fails to describe the color-suppressed 
weak decay modes. For example, the predicted decay rate of D° — > K 71° by naive factoriza- 
tion is too small by two orders of magnitude compared to experiment. Second, the hadronic 
matrix element under factorization is renormalization scale \x independent as the vector or 
axial- vector current is partially conserved. Consequently, the amplitude Cj(/x)(0)f act is not 
truly physical as the scale dependence of Wilson coefficients does not get compensation from 
the matrix elements. 

A plausible solution to the aforementioned scale problem is to extract the \x dependence 
from the matrix element (O(fi)), and combine it with the //-dependent Wilson coefficients to 
form //-independent effective Wilson coefficients. After making a physical amplitude explic- 
itly //-independent, the factorization hypothesis is applied to the hadronic matrix elements. 
However, the //-evolution factor extracted from (0(//)) depends on an infrared cutoff, which 
is originally implicit in (O (//)). Since an off-shell external quark momentum is usually chosen 
as the infrared cutoff, the //-evolution factor also contains a gauge dependent term accompa- 
nied off-shell external quarks. Therefore, this solution, though removes the scale and scheme 
dependence of a physical amplitude in the framework of the factorization hypothesis, often 
introduces the infrared cutoff and gauge dependence. 

In this paper we shall show that the above controversies can be resolved by perturbative 
QCD (PQCD) factorization theorem. In this formalism, partons, i.e., external quarks, are 
assumed to be on shell, and both ultraviolet and infrared divergences in radiative corrections 
are isolated using the dimensional regularization. Because external quarks are on shell, 
gauge invariance of the decay amplitude is maintained under radiative corrections. The 
obtained ultraviolet poles are subtracted in a renormalization scheme, while the infrared poles 



are absorbed into nonperturbative bound-state wave functions. Various large logarithms 
produced in radiative corrections are summed to all orders into the Wilson and Sudakov 
evolution factors. The remaining finite piece is grouped into a hard decay subamplitude. 
The decay rate is then factorized into the convolution of the hard subamplitude, the Wilson 
coefficient, and the Sudakov factor with the bound-state wave functions, all of which are 
well-defined and gauge invariant. The partition of the nonperturbative and perturbative 
contributions is quite arbitrary. Different partitions correspond to different factorization 
schemes. However, the decay rate, as the convolution of the above factors, is independent of 
factorization schemes as it should be. 

In Sec. II we review the conventional solutions to the scale and scheme dependence 
present in the factorization hypothesis, and their problems. Gauge invariance of radiative 
corrections is explicitly justified to all orders in Sec. III. The PQCD approach is introduced 
in Sec. IV. Explicit calculations of the evolution factor gi([i) to be defined below are shown 
in Sec. V. Section VI is the conclusion. 

II. GAUGE DEPENDENCE AND INFRARED SINGULARITY 

The aforementioned scale problem with naive factorization can be circumvented in two 
different approaches. In the first approach, one incorporates nonfactorizable effects into the 
effective coefficients [|]-||: 

af = Cl (/i) + cM (— + X i(j4) , of = C M + ci(/i) (jf + X 2 (/i)) , (3) 

where nonfactorizable terms are characterized by the parameters \i- F° r the decay B —>■ 
D + 7i~ , xi is given by 

C2 (/i)xi(/i) = (*(/*) + C 4^) 4 BD "V) + c 2 (/,)4 BA7r V) , (4) 

where 
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(BD,n) = (£> + 7T \(cb) v _ A (du) v _ A \B°) 

" (D+\(cb) v _ A \B°)(n-\(du) v _ A \0) 



£ (bd,«) = (D + n-\l(c\ a b) v „ A (d\ a u) v „ A \W) 
(D+\(cb) v _ A \B°)(K-\(du) v _ A \0) ' 

are nonfactorizable terms originated from color-singlet and color-octet currents, respectively, 
(giA a g 2 )^_ A = 9iA a 7 M (l — 75)92- The \i dependence of Wilson coefficients is assumed to be 
exactly compensated by that of Xi(//) @]. That is, the correct /1 dependence of the matrix 
elements is restored by the nonfactorized parameters Xiil 1 )- However, there are two potential 
problems with this approach. First, the renormalized 4-quark operator by itself still depends 
on /i, though the scale dependence of (0(h)) is lost in the factorization approximation. 
Second, to the next-to-leading order (NLO), the Wilson coefficients depend on the choice of 



the renormalization scheme. It is not clear if Xi(^) can restore the scheme dependence of 
the matrix element. 

In the second approach, it is postulated that (0(pb)) is related to the tree-level hadronic 
matrix element via the relation (0((j,)) = g(fj,)(0)t ree and that g(p) is independent of the 
external hadron states. Then schematically we can write 

c(jj)(0(jj)) = c(/i)g(iJ,)(0) tTee = c eff (0) tree . (6) 

The factorization approximation is applied afterwards to the hadronic matrix element of the 
operator O at the tree level. Since the tree-level matrix element (O)tree is renormalization 
scheme and scale independent, so are the effective Wilson coefficients df and the effective 
parameters af 1 expressed by ||[| 

< = of + cf (±- + Xl ) , af = cf + cf (±- + X2 ) . (7) 

However, the problem is that we do not know how to carry out first-principles calculations 
of (O(fjL)) and hence g(fjL). It is natural to ask the question: Can g(p) be calculated at the 
quark level in the same way as the Wilson coefficient c(/i) ? One of the salient features 
of the operator product expansion (OPE) is that the determination of the short-distance 
c(/i) is independent of the choice of external states. Consequently, we can choose quarks as 
external states in order to extract c(/x). For simplicity, we consider a single multiplicatively 
renormalizable 4-quark operator O (say, + or OJ) and assume massless quarks. The QCD- 
corrected weak amplitude induced by O in full theory is 
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(0) q , (8) 



where 7 is an anomalous dimension, p is an off-shell momentum of the external quark lines, 
which is introduced as an infrared cutoff, and the non-logarithmic constant term a in general 
depends on the gauge chosen for the gluon propagator. The subscript q in (Q) emphasizes 
the fact that the matrix element is evaluated between external quark states. In effective 
theory, the renormalized (0(h)) q is related to (0) g in full theory via 



<0(aO> s 




2 - P 
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+ r 
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(o)< 



= gi(^-p\\)(0) v (9) 



where g' indicates the perturbative corrections to the 4-quark operator renormalized at the 
scale /i. The constant term r is in general renormalization scheme, gauge and external 
momentum dependent, which has the general expression 0: 



,NDR,HV , x A 



+ Xr\ (10) 

where NDR and HV stand for the naive dimension regularization and 't Hooft-Veltman 
renormalization schemes, respectively, and A is a gauge parameter with A = corresponding 



to Landau gauge. Matching the effective theory with full theory, Atuii = A c s = c(fi){0(fi)) q , 
leads to 



4tt V 2 fi 2 



c{n)= 1 + ^ -Iln^ + d , (11) 



where d = a — r. Evidently, the Wilson coefficient is independent of the infrared cutoff and 
it is gauge invariant as the gauge dependence is compensated between a and r. Of course, 
c(/i) is still renormalization scheme and scale dependent. 

Since A eS in full theory [Eq. (8)] is \x and scheme independent, it is obvious that 

c' eff = c(/iy (/i, ~V\ A) (12) 

is also independent of the choice of the scheme and scale. Unfortunately, c' eS is subject to the 
ambiguities of the infrared cutoff and gauge dependence, which come along with g' extracted 
from (0(/i)) g . As stressed in JJJ, the gauge and infrared dependence always appears as long 
as the matrix elements of operators are calculated between quark states. Therefore, it is 
unreliable to define the effective Wilson coefficients by applying the existing calculations in 
the literature. The reason has been implicitly pointed out in M that "off-shell renormalized 
vertices of gauge-invariant operators are in general gauge dependent" . 

The existing problems associated with the off-shell regularization scheme are as follows: 

1. When working with off-shell fermions, there exists the so-called P operator || e.g., 
f(l — 75) ®^(1 — 75) which cannot be removed by the equation of motion. 

2. The finite terms are external momentum dependent (see Fig. 3 of ||) and they are 
obtained in some specific condition. For example, two incoming fermion legs 1,2 and two 
outgoing legs 3,4 with external momentum p are chosen in Figs. 3a and 3c of ||[10|, while 
legs 1,3 incoming and 2,4 outgoing with p in Fig. 3b. 

Hence we cannot avoid the gauge problems if adopting off-shell fermions and the finite 
parts of c' eS are not well-defined. To circumvent this difficulty, we should work in a physical 
on-shell scheme and employ the dimensional regularization for infrared divergences. Gauge 
invariance of the decay amplitude is maintained under radiative corrections, and the infrared 
poles are absorbed into the hadronic matrix element as stated in the Introduction. Conse- 
quently, the effective coefficient c cS = c(p)g(p) does not suffer from the gauge ambiguity. 

III. GAUGE INVARIANCE IN ON-SHELL REGULARIZATION 

In this section we show that gluon exchanges among the on-shell quarks involved in heavy 
meson decays, including the spectator quarks, indeed give gauge invariant contributions. We 
present the proof in the covariant gauge d ■ A = 0, in which the gluon propagator is given 
by (-i/l 2 )N^(l) with 

7V^(0 = <r-(l-A)^, (13) 
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FIG. 1. Graphical representation for the proof of Eq. (14). 



where A is the gauge parameter. We shall show that the quark amplitude A{ u n with the 
spectator quarks included are independent of A to all orders, namely, 



A 



dA 



full 



dX 







The differential operator applies only to gluon propagators, leading to 

IH V 



dX 



X- 



P 



v a [FN ca, + V , N afi ] 



(14) 



(15) 



with the special vertex v a = l a /(2l 2 ). The loop momentum Z M (l v ) carried by the differen- 
tiated gluon contracts with a vertex in A full , which is then replaced by the special vertex 
v a . Eq. flT5|) is graphically described by the first expression in Fig. 1, where the arrrow 
represents Z M {l u ) contracting with the gluon vertex, and the square represents v a . 

The contraction of Z M (l u ) leads to the Ward identity jTIJ shown in the second expression 
of Fig. 1, where the solid lines may represent quarks or gluons. Summing all the diagrams 
with various differentiated gluons, those embedding the special vertices cancel by pairs. For 
example, the pair cancellation occurs between the first and last diagrams in the second 
expression of Fig. 1. Only the diagram, in which the special vertex moves to the outer end of 
the quark line, is left. This diagram comes from the second term in the following expression, 



i( / k+ /+M) 
(k + I) 2 - M 2 



—i f)u(k) = u(k) 



fi+ f+M 
(k + I) 2 - M 2 



{fc-M)u{k) 



(16) 



where u is the fermion spinor associated with an external quark. The first term is canceled 
by the term from the contraction of I with the adjacent vertex. If all the external quarks 
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FIG. 2. Vertex corrections to the 4-quark operators 0\ and 2 . 

are on shell, the second term vanishes because of the equation of motion {Jk — M)u(k) = 0. 
Then we arrive at the desired result (|T4|). 

We take one-loop corrections as an example to elucidate the above proof. We consider 
only the gauge-dependent part of the gluon propagator [see Eq. ([131) 1, 



— i 



A) 



IH 1 



f- 



(17) 



in the loop calculations and demonstrate that the result vanishes after summing all the 
diagrams. The gauge-dependent part of Fig. 2(a) reads 



.G f 



/r 9e = *;^cKM(]--- A ),/:/, 



2,.e 



X 



dPk u 3 t a ^(l - 75 ) ( A+ fr + M x ) fax u^Yil ~ 7 5 )(i>2- h + M 2 ) fa 2 



(27T) 



D 



k%k + Pl y - Mi)[{k - P2 y - Mi) 



where Vckm is the relevant Cabibbo-Kobayashi-Maskawa matrix elements. 
To proceed, we replace the Jk which is adjacent to Ux as 

h =A+ fr - Mx - (A - Mx) , 

and the Jk adjacent to u 2 as 

fi = - fa+ fi + M 2 + (jf> 2 -M 2 ). 
Applying Eq. (|l^) and the equation of motion, (/> — M)u = 0, Eq. (|Tj|) becomes 



(19) 



(20) 



4 =^ckm( v 1-aj^ y ^ jfci • ( 21 ) 



Likewise, one can apply the same trick to the calculations of Figs. 2(b)-(f) and obtain 

jgauge _ _jgauge _ _j90,uge __ jgauge J 



jgauge jgauge 

G F .. / \ o ,„ f d D k 1 



/ — Vckm( 1 - X)g 2 s H e C F J T^vd^Z^^ 1 - 7s)ui ^47 M (1 - 75)^2 



= -^Vckm(i - A )^(t- " 7")«37m(1 -75)«i «47"(1 -75)«2, (22) 

where e uv(IR) = 4 — .D is the ultraviolet (infrared) pole, and Cf = (TV,? — 1)/(2JV C ) with JV C be- 
ing the number of colors. After lengthy but straightforward calculation, the renormalization 
constant of a fermion with mass M is found to be 

z 2 = 1 — — G,f 37s + 4 — 3 In 



4tt Ve uv ' ' ' 4vr 2 /i 2 e IR . 

+(1 _ a) ^ Cf (A_A). (23) 

We see that, contrary to the gauge-independent part of Z 2 — 1, the gauge-dependent contri- 
bution due to the fermion wavefunction renormalization 



V2 



< ckm( 1 - X )^ C p(— ~ — )«37a.(1-75)«i «47 m (1 -7s)«2 (24) 



'UV ^IR " 



is free of mass singularity. Summing over all the contributions, it is obvious that the final 
result indeed vanishes as it should be. 

IV. PQCD FACTORIZATION THEOREM 

We have shown that radiative corrections to a decay amplitude of on-shell external quarks 
are gauge invariant to all orders. The one-loop diagrams have been evaluated explicitly, 
whose results confirm our proof. Next we shall explain how to treat the infrared poles in the 
PQCD factorization theorem. The one-loop contributions in full theory are ultraviolet finite 
because of the current conservation stated above. The existence of the infrared poles simply 
signifies the nonperturbative dynamics, which demands the inclusion of bound-state effects 
into the formalism of heavy meson decays. The standard treatment of infrared poles is to 
absorb them into a universal meson wave function. To absorb the infrared poles associated 
with the b quark, such as those from the self-energy corrections, it is necessary to introduce 
a B meson wave function. That is, we must take into account the spectator quark of the B 
meson in order to develop a complete theory of heavy meson decays. 

Accordingly, the decay amplitude Af u u to one-loop in full theory can be rewritten as 



A ^^ as c , i Mw , 'i M '/ 
Ami = 1 + — h 7 In -rj- + 7 In ha 

An \e IR M b n f j 

a s f M w M b N 



+0(a 

47re IR y 



2 ) 



(25) 



where the factorization scale /u,f arises from the dimensional regularization of infrared diver- 
gences, and the factorization of the infrared pole is performed in the minimal subtraction 
scheme. The anomalous dimensions of the logarithms ln(Mw/M b ) and ln(M&//4/), 7 and 
7', respectively, are different, since the latter involves an extra contribution related to the 
spectator quark. The factor containing the infrared pole can be formulated as a matrix ele- 
ment of a nonlocal operator, which is the definition of a meson wave function 0(/i/). A wave 
function, describing the amplitude that a parton carries a fraction of the meson momentum, 
cannot be derived in perturbation theory. It must be parametrized as a function of parton 
momentum fraction. 

We further factorize the infrared finite part into 
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full 
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1 + 47l v 7ln A4 +7l V + a_a , 



(26) 



The first factor, characterized by the matching scale Mw, is identified as the Wilson coeffi- 
cient c((jl) after summing ln(M H // / u) to all orders using renormalization group (RG) equations. 
The second factor, characterized by the b quark mass M b , is the hard subamplitude which 
will be denoted by H(M b , fj,, fj,f) below. Extending the above procedures to all orders, we 
obtain the factorization formula for B meson (not b quark) decays, 



A 



full 



c(fjL)H(M b ,fi,fi f )^(jjL f ), (27) 

which is exactly the three-scale factorization formula for exclusive nonleptonic decays derived 



in [[TI|. Note that Eq. (p7|) in fact denotes a convolution relation, because the momentum 
fractions should be integrated out. 

Compared to Eq. (H), the matrix element (0(h)) corresponds to 



(0(fi)) = H(M b ^,fi f )^ f ) 



(28) 



Summing ln(fi/M b ) in H to all orders using RG equations, we obtain an evolution factor 
gi(fj), whose behavior from \x to M b is governed by the same anomalous dimension as that 
of c(\i). Summing \n(M b /(j,f) in H to all orders, we obtain another factor g 2 (fif) describing 
the evolution from M b to /i/, whose anomalous dimension differs from that of c(/x) because 
of the inclusion of the dynamics associated with spectator quarks. Note that g 2 is part of 



the Sudakov evolution obtained in |[L2| . Hence, the /1 dependence of H is extracted as 



H(M b , /1, fi f ) = gi (n)g 2 (n f )H(M b , M b , M b 



(29) 



The combination of c, gi, and g-i leads to the effective coefficient 



c eff = c(fi) gi (fi)g 2 (fx f ) , (30) 

which is not only \x and scheme independent but also gauge invariant. The factor g{p) in 
Eq. (|B|) can be identified as gi(/i)g 2 (/!/), which describes the evolution down to the factor- 
ization scale. However, g\{p)gi{p,f) contains a matching condition at the scale M& between 
the Wilson and Sudakov evolutions with different anomalous dimensions. Therefore, there 
is an ambiguity of the matching condition: the two evolutions can also match at rM b with 
r a constant of order unity. Obviously, Eq. ( |30f) is subtler than the naive definition of c in 
Eq. (§. 

The matrix element (0) tre e m Eq. (|) is then identified as 

(O)tree = H trcc (M bl M h M b )<f>(jl f ) , (31) 

where the hard subamplitude is evaluated to lowest order with one hard gluon exchange, 
since all large logarithms have been organized by RG equations. We emphasize that the 
factorization hypothesis for (0) tre e i n the conventional approach is not necessary in the 
PQCD formalism. The purpose of the factorization hypothesis is to simplify the decay 
amplitude into products of decay constants and form factors, which are then parametrized 
as various models. To have a better fit to experimental data, nonfactorizable contributions, 
parametrized as x [ see Eq. (7)], are included. Note that i^tree in the PQCD approach 
includes both factorizable contributions (form factors), when the hard gluon attaches to the 
two quarks in a meson, and nonfactorizable contributions (octet amplitudes), when the hard 
gluon attaches to the quarks in different mesons. Therefore, we may compute all possible 
diagrams for i^trce HH and convolute them with the same meson wave functions <fi. That is, 
we use the single parametrization, i.e., the meson wave functions, for both factorizable and 
nonfactorizable contributions based on Eq. fl3"T|). In this sense the PQCD formalism is more 
systematic. 

At last, we explain how to handle the nonperturbative meson wave functions with the 
dependence of the factorization scale fif. It can be shown that these wave functions are 
universal for all decay processes involving the same mesons. For example, the B meson wave 
function for the nonleptonic decays B — > D^ir(p) and for the radiative decay B — ► K*^ 
is the same. This universality can be easily understood, since a wave function collects 
long-distance (infrared) dynamics, which should be insensitive to short- distance dynamics 
involved in the decay of the b quark into light quarks with large energy release. Based on 
the universality of wave functions, the application of factorization formulas is as follows 
r3fl . We evaluate the Wilson and Sudakov evolutions down to a factorization scale /i/ and 



the hard subamplitude for a decay mode, say, B — > K*j, in perturbation theory. These 
calculations are simply performed at the quark level with infrared poles dropped (in the 
minimal subtraction scheme). Adjust the B meson wave function such that the predictions 
from the relevant factorization formula match the experimental data. At this stage, we 
determine the B meson wave function defined at the scale /i/. Then evaluate the Wilson and 
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Sudakov evolutions down to the same scale /z/ and the hard subamplitude for another decay, 
say, B — > Dix . Convolute them with the same B meson wave function and make predictions. 
At this stage, there are no free parameters in the formalism. With the above strategy, the 
PQCD factorization theorem possesses predictive power. 

The main uncertainties in the PQCD factorization theorem come from higher-order cor- 
rections to the hard subamplitude and higher-twist corrections from the Fock states other 
than the leading one with only valence quarks which we are considering here. Accord- 
ing to Eq. fl29|), the argument of the running coupling constant in the hard subamplitude 
H(M b , M b , M b ) should be set to the b quark mass M b , implying that the next-to-leading- 
order diagrams give about a s (M b )/it ~ 10% corrections. Since H is characterized by M b , 
the next-to-leading-twist correction from the Fock state with one more parton entering the 
hard subamplitude is about ftf/M b ~ 10%. Note that meson wave functions are usually 
defined at the factorization scale fif ~ 0.5 GeV [p!4| . We believe that other nonperturbative 
corrections, such as final-state interactions, should play a minor role because of the large 
energy release involved in two-body B meson decays. If the hadronic matrix elements are 
evaluated using the factorization approximation (i.e. vacuum insertion approximation), the 
related uncertainties have been discussed in length in |T7]. 



In conclusion, all the factors in the PQCD formalism are well-defined (including the 
nonperturbative meson wave functions) and gauge invariant. Physical quantities obtained in 
this formalism are scale and scheme independent. We have applied this approach to exclusive 
semileptonic, nonleptonic, and radiative B meson decays and the results are very successful. 
Nonfactorizable contributions have been calculated, and found to play an important role in 
the decays B — > J/ipK^ [P~5| . The opposite signs of a 2 /a± in bottom and charm decays 
have been explained by the effects of the Wilson evolution fl5| . The mechanism for the 
sign change of the nonfactorizable contributions in bottom and charm decays have also been 
explored [[U| , which is closely related to the success and failure of the large- A c limit in charm 
and bottom decays, respectively. 

V. EFFECTIVE WILSON COEFFICIENTS 

In this section we present the results for the evolution factor g\ (/i) which describes the 
evolution from the scale /i to M b for the current-current operators 0\ and 2 -* Setting 
;Uj = M b , the effective Wilson coefficients obtained from the one-loop vertex diagrams Figs. 
2(a)-(f) for the operators O; have the form: 



„eff 



M =cM + ^U^\n^ + A aM, 

H= M h 4:71 \ fl J u 



*The complete results of gi(/J,) for AB = 1 transition current-current operators 0±,02, QCD- 
penguin operators O3, • • • , Oq and electroweak penguin operators Oj, • ■ ■ , 0\q are given in |TT 
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.off 



C, 



M =cM + ^( 1 ^ln^+rA c^), (32) 



where the superscript T denotes a transpose of the matrix, and the anomalous dimension 
matrix 7^ due to the one-loop vertex corrections has the well-known expression 

+-(? _y. (33, 

The matrix r gives momentum-independent constant terms which depend on the treatment of 
75. Working in the (massless) on-shell scheme and assuming zero momentum transfer squared 
between color-singlet currents, i.e. (pi — P3) 2 = as well as (pi + P2) 2 = (— P2 + P3) 2 ~ rn 2 
for 0\ operators and (pi — p 4 ) 2 = 0, (pi + p 2 ) 2 = (— P2 + Pi) 2 ~ Tn 2 for 2 operators (see 
Fig. |2| for momentum notation), we obtain 

rNDR= (-9 _ 3 9 )' rHV= (-7 ~l) (34) 

in NDR and HV schemes, respectively. It should be accentuated that, contrary to the 
previous work |5[] 

4Br = (^_ 3 7 z), 4v° = y_ 5 7 ), (35) 

obtained in Landau gauge and off-shell regularization, the matrix r given in (133) is gauge 
invariant ! 

Two remarks are in order. First, there are infrared double poles, i.e., 1/e 2 , in the 
amplitudes of Figs. 2(a)-2(d), but they are canceled out when adding all amplitudes together. 
Second, care must be taken when applying the projection method to reduce the tensor 
products of Dirac matrices to the form r (g> T with T = 7 M (1 — 75). For example, a direct 
evaluation of the tensor product 7 a faiT <g> Y p 2 l a yields (e = 4 — D) 

7a far ® r M a = -c(pi • P2)T ® r (36) 

in the NDR scheme with the on-shell condition being applied first to the massless quarks 
followed by Fierz transformation, whereas the projection method of |T^JT0fl] leads to 



la far ® r fr 2l a = ^ • p 2 )r <g> r + e, (37) 

where E stands for the evanescent operator (EO). This means that it is incorrect to take 
the coefficient of T <g> T in Eq. (^) directly without taking into account the effect of EOs. 
Note that we have applied Eq. ([3(]) to show the absence of infrared double poles in the total 
amplitude. 

In order to check the scheme and scale independence of c° s , it is convenient to work in the 
diagonal basis in which the operators 0± = \{0\ ± 2 ) do not mix under renormalization. 
Then (see e.g. [^] for the general expression of c(ix)) 
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(38) 



where c± = ci ± C2, /3o — 11 — gTi/ with n/ being the number of flavors between My/ and 
pi scales, B± specifies the initial condition of c{m w ): c{m w ) = 1 + ^ i?± and it is 
75-scheme dependent, and J± = 7i° ) /V( 2 / ? o) _ fy^/(2/3o), with f3 1 = 102 - 38n//3. The 
scheme-dependent anomalous dimensions j± are given by p0|,|7|: 
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(39) 



where 7± are the two- loop anomalous dimensions of O-t, 7j is the anomalous dimension of 
the weak current in full theory, and the parameter k± distinguishes various renormalization 
schemes: K± = in the NDR scheme and K± = =F4 in the HV scheme. As shown in 
|p0| , B± — J± is 75-scheme independent. Therefore, the effective Wilson coefficients c^f are 
scheme independent if we are able to show that (r T + J) is independent of the choice of 
renormalization scheme. Since the short-distance Wilson coefficients are independent of the 
choice of external states, one can show the independence of 7^ — 27J from external states 
|20fl . In the on-shell scheme, jj vanishes up to the two- loop level. It follows from Eq. ([34]) 
that r + = —6, —7 + 7/3 and r_ = 12, 7 + 7/3 in NDR and HV schemes, respectively, with 
fermions being on-shell. Then it is easily seen that r T + J is indeed renormalization scheme 
independent. To the leading logarithmic approximation, 
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Hence, the scale independence of the effective Wilson coefficients follows. 

Since the weak current is partially conserved, its anomalous dimension 7j is zero. How- 
ever, if fermions are off-shell, 7j is non-vanishing at the two-loop level in the HV scheme. 
To maintain the requirement that 7j = 0, one can force a vanishing jj in this case by 
applying a finite renormalization term to the weak current. (Note that in this new choice, 
7± = 7± , B'± = B± — 7j//3o, and B'± — J'± is still scheme independent.) Using the identity 
7^ (on-shell) =7 ( - 1 - ) (off-shell) —2'jj, we find that 7^ 1 - ) in the off-shell fermion scheme is given 
bypj: 
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From Eqs. fl55|) and (PD, it is straightforward to show that r± + J± is 75-scheme indepen- 
dent in the off-shell regularization. As a result, c eff is renormalization scheme independent, 
irrespective of the fermion state, on-shell or off-shell. 

VI. CONCLUSION 

In this paper we have shown how to construct a gauge invariant and infrared finite theory 
of exclusive nonleptonic B meson decays based on PQCD factorization theorem. Gauge in- 
variance is maintained under radiative corrections by working in the physical on-shell scheme. 
The infrared divergences in radiative corrections should be then isolated using the dimen- 
sional regularization. The resultant infrared poles are absorbed into the universal meson 
wave functions, which can be determined once for all from experimental data. The absorp- 
tion of the poles associated with the b quark requires the inclusion of the spectator quark 
into the theory. The remaining finite contributions form a hard subamplitude. Applying RG 
analyses to sum various large logarithms in the above factorization formula, the scale and 
scheme dependences are removed. Hence, in the PQCD formalism physical quantities are 
guaranteed to be gauge invariant, infrared finite, scale and scheme independent. By work- 
ing out the evolution factor gi(fi) explicitly, we have constructed gauge invariant, scale and 
scheme independent effective Wilson coefficients c(/i)gi(/i) at the factorization scale /// = M},. 
We have shown explicitly that c eff are renormalization scheme and scale independent. 

We shall take one of the exclusive nonleptonic B meson decay modes as an example to 
demonstrate how to construct a factorization formula explicitly. This work will be published 
elsewhere. 
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